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Abstract 

A method for real-time estimation of parameters 
in a linear dynamic state space model was developed 
and studied. The application is aircraft dynamic model 
parameter estimation from measured data in flight for 
indirect adaptive or reconfigurable control. Equation 
error in the frequency domain was used with a recursive 
Fourier transform for the real-time data analysis. 

Linear and nonlinear simulation examples and flight 
test data from the F-18 High Alpha Research Vehicle 
(HARV) were used to demonstrate that the technique 
produces accurate model parameter estimates with 
appropriate error bounds. Parameter estimates 
converged in less than 1 cycle of the dominant dynamic 
mode natural frequencies, using control surface inputs 
measured in flight dining ordinary piloted maneuvers. 
The real-time parameter estimation method has low 
computational requirements, and could be implemented 
aboard an aircraft in real time. 
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Nomenclature 

system matrices 
expectation operator 
acceleration due to gravity, ft/sec 2 
altitude, ft 

imaginary number = -J-l 
Mach number 
total number of samples 
real part 

measurement noise covariance matrix 
time 
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V , true airspeed, ft/sec 

x, u, y state, input, and output vectors 

z, measured output vector at time i At 

a angle of attack, rad 

[1 sideslip angle, rad 

8 a , S r aileron, rudder deflections, rad 

S e , S s elevator, stabilator deflections, rad 

£>,, Kronecker delta 

y 

Vi discrete measurement noise vector 

G~ variance 

co angular frequency, rad/sec 

9 p-dimensional parameter vector 

superscripts 
T 

t 

-1 

subscripts 

i value at time i At 

o trim or initial value 


transpose 

complex conjugate transpose 
discrete Fourier transform 
estimate 
matrix inverse 


Introduction 

Real-time identification of dynamic models is a 
requirement for indirect adaptive or reconfigurable 
control 1 . One approach for satisfying this requirement 
is to assume the dynamic model has a linear structure 
with time-varying parameters to account for changes in 
the flight condition, stores, configuration, remaining 
fuel, or from various types of failures, wear, or damage. 
The task is then to identify accurate linear model 
parameter estimates from measured data in real time, so 
that the adaptive control logic can make the necessary 
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changes to the control law to achieve stability and 
performance goals. 

Two main problems plague accurate real-time 
parameter estimation: noise and data information 
content. It is difficult to design a parameter estimation 
technique that is insensitive to noise but still responds 
rapidly to sudden changes in the system dynamics, 
mainly because it takes a fairly long data record to 
distinguish noise from a sudden change in the 
dynamics. This problem can be handled in the time 
domain using recursive least squares and a "forgetting 
factor" 2 * , or by using sequential batch least squares with 
short data records and including various constraints in 
the parameter estimation cost function 3,4 . If an 
extended Kalman filtering approach 5 ' 8 is used, 
discriminating signal from noise is implemented 
through weighting matrices that represent assumed 
measurement and process noise covariances. For all of 
these time domain methods, some adjustment of one or 
more tuning parameters must be done in simulation. In 
addition, the standard errors for the model parameter 
estimates, which are important both for failure 
detection and adaptive or reconfigurable control, 
cannot be accurately and reliably computed using 
recursive or sequential batch time domain methods. 

In the context of airplane flight, lack of 
information content in the data can be problematic 
because there are frequently extended periods where 
the control and state variables are fairly constant. 

Signal levels are at or below the (relatively constant) 
noise level. In this circumstance, a time domain 
regression method will give very inaccurate parameter 
estimates unless the estimation is regularized by 
including a term in the cost function that penalizes 
movement of the parameters away from some a priori 
known values (e.g., values from wind tunnel tests), 
and/or a term that penalizes time variation of the 
parameter estimates 4 . Tuning parameters are required 
for this approach, because the magnitude of the penalty 
term(s) must be balanced properly relative to the least 
squares part of the cost function used for parameter 
estimation based on measured data. Another problem 
that falls in the category of poor data information 
content is data collinearity due to the control system 9 . 
Many control laws move more than one control surface 
at the same time, or move control surfaces in 
proportion to state variables with a small time delay. 
When states and controls are nearly proportional to one 
another, it is impossible to identify individual stability 
and control derivatives from the measured data alone. 

There are many parameter estimation methods, 
but the requirement of being simple enough to be 


implemented in real time aboard the aircraft narrows 
the field. In particular, any method that iterates through 
the data must be eliminated. The current work is an 
investigation of a single-step frequency domain method 
for the real-time parameter estimation task, and an 
evaluation of its suitability for aircraft problems. This 
real-time parameter estimation method was first 
proposed as a component of a technique for in-flight 
system identification 10 . In the present work, the 
real-time parameter estimation method is developed 
further and applied to realistic simulation and flight test 
data. 

The next section gives the problem statement and 
outlines the necessary theory. Following this, the 
real-time parameter estimation method is applied to a 
simulation example, where a linear truth model is used 
with outputs corrupted by noise similar to that seen in 
flight, and noise that is worse than usual. The 
application is identifying an accurate model for the 
longitudinal rigid body dynamics of a conventional 
fighter. The real-time parameter estimation procedure 
is further demonstrated with a longitudinal example 
using flight test data, and a nonlinear simulation 
example where both longitudinal and lateral/directional 
parameters are estimated. 

Theoretical Development 

Airplane dynamics can be described by the 
following linear model equations 11 : 

i(f) = Ax[t) + Bu(t) (1) 

*(°) = x„ (2) 

y(t) = Cx(t)+Du(t) (3) 

Z i =y i +V i i = 1,2, ...,N (4) 

Matrices A, B , C, and D in Eqs. (1) and (3) 
contain stability and control derivatives, which are 
assumed as constant model parameters to be estimated 
from flight data. Repeating the parameter estimation at 
short intervals produces piecewise constant estimates 
for time-varying model parameters in the linear model 
structure. The input quantities are control surface 
deflections (S e or S s , S a , S A, with the states selected 
from air data ( V, , a, (3 ) , body axis angular velocities 
( p,q,r ) , and Euler angles ((f), 0, if/). Output 
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quantities can include the states and translational 
accelerations [a x ,a y ,a z 


h = -z X | j°>n \ («) - A k x{n)~ B k ii (n) 

^ n — 1 


( 11 ) 


Equation Error in the Frequency Domain 

The finite Fourier transform of a signal x(t ) is 
defined by 


x{co) = J Q x(t)e > at dt (5) 

which can be approximated by 


N - 1 

x(co) ~ At ^ x k e - - 7 ®' 1 (6) 

!= o 

Subscript i indicates the variable value at time i At, and 
At is the sampling interval. The summation in Eq. (6) 
is defined as the discrete Fourier transform, 

X{a>)J£x,e-^ (7) 

i=0 


so that 


x(a>) ~ X(a>)At (8) 

Some fairly straightforward corrections 12 can be 
made to Eq. (8) to remove the inaccuracy resulting 
from the fact that Eq. (8) is a simple Euler 
approximation to the finite Fourier transform of Eq. (5). 
However, if the sampling rate is much higher than the 
frequencies of interest (as is true in this case), then the 
corrections are small and can be safely ignored. 

Applying the Fourier transform to Eqs. (1) and (3) 

gives 


jcox^co) = Ax[co) + Bu[co) 

(9) 

y(co) = C x(co) + Du(co) 

(10) 


When the states, outputs, and inputs are 
measured, individual state or output equations from 
vector Eqs. (9) or (10) can be used in an equation error 
formulation to estimate the stability and control 
derivatives contained in matrices A, B, C, and D. 


where A k and B k are the At h rows of matrices A and 
B , respectively, and x k (n) is the Ath element of vector 
x for frequency co n . Symbols x(n) and u(n ) denote 

the Fourier transform of the state and control vectors 
for frequency co n . There are m terms in the summation, 
corresponding to m frequencies of interest, and each 
transformed variable depends on frequency. Similar 
cost expressions can be written for individual output 
equations from vector Eq. (10). Denoting the vector of 
unknown model parameters in A k and B k by 9 , the 
problem can be formulated as a standard least squares 
regression problem with complex data. 


Y = X9 + £ 


where 


Y = 


je->\ wO) 

j°> 2 M 2 ) 


j(O m X k (m) 


i r (l) 

5 r (l) 

x T (2) 

u r (2) 

x T (m ) 

u T (m) 


( 12 ) 


(13) 


(14) 


and £ represents the complex equation error in the 
frequency domain. The least squares cost function is 

J = ^(Y-X9) t (Y-X9) (15) 

which is identical to the cost in Eq. (11). The 
parameter vector estimate that minimizes this cost 
function is computed from 13 

() ’av(.Y .Y)' AV(.V }’) (16) 

The estimated parameter covariance matrix is 


For the At h state equation of vector Eq. (9), the 
cost function is 


cov(6>) = E{{e-e){e-e) r | = o 2 \Re{x f z)] '(17) 
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where the equation error variance a 1 can be estimated 
from the residuals. 


a 1 - 


(m-p) 


{y -xe) f {y -xe) 


(18) 


and p is the number of elements in parameter vector 6 . 
Parameter standard errors are computed as the square 
root of the diagonal elements of the cov{0 j matrix from 
Eq. (17), using a 2 from Eq. (18). 


Recursive Fourier Transform 

For a given frequency co, the discrete Fourier 
transform in Eq. (7) at sample time i is related to the 
discrete Fourier transform at time i — 1 by 


0.04 Hz on the interval [0.1 -1.5] Hz was found to be 
adequate, giving 36 evenly spaced frequencies for each 
transformed time domain signal. Excluding zero 
frequency removes trim values and measurement 
biases, so it is not necessary to estimate bias 
parameters. Using a limited frequency band for the 
Fourier transformation confines the data analysis to the 
frequency band where the system dynamics reside, and 
automatically filters wide band measurement noise or 
structural response outside the frequency band of 
interest. 

For airplane dynamic modeling, the number of 
time domain signals to be transformed is usually low (9 
or less - more if there are many control surfaces), so 
this approach requires a small amount of computer 
memory. Since the data analysis is done in the 
frequency domain, the memory required is fixed and 
independent of the time length of the flight maneuvers. 


Xi((o) = + x t e~ jwi ^ (19) 


where 


e ~jcoiAi _ e -jwAl e -j<o(i-l)At (20) 

The quantity e - - 7 ® 4 * ' is constant for a given 
frequency and constant sampling interval. It follows 
that the discrete Fourier transform can be computed for 
a given frequency at each time step using one addition 
in Eq. (19) and two multiplications - one in Eq. (20) 
using the stored constant for frequency co, and 

one in Eq. (19). There is no need to store the time 
domain data in memory when computing the discrete 
Fourier transform in this way, because each sampled 
data point is processed immediately. Time domain data 
from all preceding maneuvers can be used in all 
subsequent analysis by simply continuing the recursive 
calculation of the Fourier transform. In this sense, the 
recursive Fourier transform acts as memory for the 
information in the data. More data from more 
maneuvers improves the quality of the data in the 
frequency domain without increasing memory 
requirements to store it. In addition, the Fourier 
transform is available at any time i At . The 
approximation to the finite Fourier transform is 
completed using Eq. (8). 

Rigid body dynamics of piloted aircraft lie in the 
rather narrow frequency band of approximately 
0.01-1.5 Hz. It is therefore possible to select 
closely-spaced fixed frequencies for the Fourier 
transform and the subsequent data analysis. For all the 
examples studied in this work, frequency spacing of 


The recursive Fourier transform update need not 
be done for every sampled time point. Skipping some 
time points effectively decimates the data prior to 
Fourier transformation. This saves computation, and 
does not adversely impact the frequency domain data 
because the Nyquist frequency (equal to Vi the 
sampling frequency) is usually much higher than the 
relatively low frequencies being used in the recursive 
Fourier transform. 


Examples 

For longitudinal aircraft short period dynamics, 
the state vector x and input vector u in Eq. (1) are 
defined by 


x = [a q] T M = [c> e ] (21) 


System matrices containing the model parameters are: 


A = 


Z a 

M a 




( 22 ) 


The above model assumes a effects can be subsumed 
into the Z' and derivatives. Parameter Z' includes 

the inertial term, i.e., Z ' = I + Z . . I n this and all the 

Lj q 

other examples, state equations were used for the 
equation error parameter estimation. 

In the first example, a perturbation elevator input 
was applied to a known linear model to produce 
simulated state and output responses. Figure 1 shows 
the elevator input 8 e . The first 8 seconds of the 
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elevator input were taken from measured flight test data 
for a longitudinal tracking task. Constant elevator 
deflection is held for the final 7 seconds, to simulate 
low data information content. The simulated aircraft is 
a conventional F-16 11 ’ 14 with forward c.g. position 
(0.2 c) in straight and level trimmed flight at 10,000 ft, 
7 deg angle of attack, and Mach 0.37. The simulated 
outputs were corrupted with 20% Gaussian random 
white noise. This made the signal-to-noise ratio 5-to-l 
for each simulated output measurement. Figure 2 
shows the simulated measured perturbation angle of 
attack a and pitch rate q. The elevator input was 
assumed to be measured without noise, which is a close 
approximation to reality. Model parameter values used 
to generate the simulated test data, called true values, 
are given in column 2 of Table 1. Parameter estimation 
was done in real time using equation error in the 
frequency domain applied to the two state equations, 
with the Fourier transform computed recursively, as 
described previously. All angular quantities were 
expressed in radians for the data analysis, but were 
plotted in degrees. Angular rates were expressed in 
radians per second for the data analysis, and plotted in 
degrees per second. 

Figure 3 shows a time history of M a model 
parameter estimates based on the simulated noisy data 
only. Plots for the other model parameters were 
similar. The model parameter estimates and standard 
errors were computed once a second using Eqs. 
(16)-(18) with frequency domain data from the 
recursive Fourier transform in Eqs. (19) and (20). This 
update rate for the parameter estimation was used 
throughout the examples, but was chosen arbitrarily. 
Parameter estimation updates can be done at a faster or 
slower rate, with the upper limit defined by the rate 
used in the recursive Fourier transform. The algorithm 
required no starting values for the parameters, and the 
estimates were not regularized in any way with a priori 
values or constraints on temporal changes in the 
parameter estimates. The first parameter estimates are 
shown at 2 seconds, because the parameter estimates 
after 1 second were poor with very large standard 
errors, due to the lack of information content in the data 
during the first second of the maneuver. Figure 3 
shows that the parameter estimates converge to the true 
value. The calculated standard errors are representative 
of the estimated parameter accuracy throughout the 
maneuver, and do not suffer from the covariance 
wind-up problem characteristic of recursive time 
domain methods. Column 3 of Table 1 contains 
parameter estimates with corresponding standard errors 
in parentheses below. These results are for the end of 
the 15 second maneuver. Every parameter estimate is 


within +1 standard error of the true value, indicating 
that the parameter estimation is accurate, and the 
estimated standard errors properly represent the true 
accuracy of the parameter estimates. The calculated 
standard errors accurately convey information on the 
quality of the parameter estimates throughout the 
maneuver, and do not become smaller with increasing 
maneuver time when there is no information in the data. 

Figure 4 shows the same simulated model pitch 
rate output using the same input to the same model, but 
with the added Gaussian random noise level raised 
from 20% to 50%, reducing the signal-to-noise ratio to 
2-to-l. In addition, two simulated data dropouts with 
values of -100 were added to the simulated pitch rate 
output. Plots in Figure 5 show that the parameter 
estimates in this case again converge to the true values, 
although the standard error values are generally higher, 
due to the increased noise level. Plots for the model 
parameters not shown were similar. The convergence 
rate of the parameter estimates to the true values was 
similar to the lower noise case, requiring approximately 
4 seconds of data. This corresponds to approximately 
1 .4 periods of the short period natural frequency for the 
simulation model. Considering that no substantial 
information is contained in the data for the first two 
seconds, this is an excellent result. Discounting the 
first two seconds, the parameters were accurately 
estimated from approximately 2 seconds of data, 
corresponding to 0.7 cycles of the short period mode. 
The real-time parameter estimation algorithm is robust 
to measurement noise levels and infrequent data 
dropouts because of the automatic filtering inherent in 
using a limited bandwidth for the recursive Fourier 
transform. In effect, the data dropouts look the same as 
high frequency noise. Column 4 of Table 1 gives the 
parameter estimates and standard errors for the 50% 
noise case at the end of the 15 second maneuver. As 
before, every parameter estimate is within +1 standard 
error of the true value, indicating that the parameter 
estimates and standard errors are accurate. 

The linear simulation and the real-time data 
analysis were programmed and run in MATLAB 5.3 15 . 
Sampling rate for the data was 40 Hz, and the recursive 
Fourier transform updates were done at 20 Hz. The 
real-time estimation algorithm ran roughly 10 times 
faster than real time (1.5 seconds for a 15 second 
maneuver) on a Gateway 450 MHz E-4200 serial 
processor, running Microsoft Windows NT 4.0. 

The next example used real flight test data from 
the F-18 High Alpha Research Vehicle (HARV) to 
demonstrate the real-time parameter estimation method. 
Figure 6 shows the measured stabilator deflection for 
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this 14 second maneuver. Measured outputs are shown 
in Figure 7. The maneuver was a sequence of doublets 
implemented by the pilot from a steady trim condition 
at 20 deg angle of attack, 24,100 ft altitude, and Mach 
0.37. Parameter estimates are plotted as a function of 
time in Figure 8, along with dashed lines indicating the 
±2<7 error bounds. Parameter estimates from standard 
batch time domain output error parameter estimation 16 
are plotted as solid lines for comparison. The real-time 
parameter estimation algorithm produced parameter 
estimates and standard errors that were in agreement 
with the batch time domain estimates after about 5 
seconds, of which 1 .5 seconds was steady trim with no 
information in the data. One cycle of the short period 
mode using the batch time domain parameter estimates 
was 7.5 seconds. Table 2 contains results at the end of 
the 14 second maneuver from batch time domain and 
real-time frequency domain parameter estimation. 

As in the simulated data cases, the standard errors 
for the parameter estimates computed by the real-time 
parameter estimation algorithm were consistent with the 
accuracy of the parameter estimates throughout the 
maneuver. Specifically, the error bounds were large at 
the beginning of the maneuver, and smaller as more 
information was obtained from the data. Except for a 
few instances that are to be expected because of 
statistical variation, the standard errors were 
representative of the accuracy of the estimated 
parameters. 

Sampling rate for the flight test data was 50 Hz, 
and the recursive Fourier transform updates were done 
at 25 Hz. The same implementation and computer as 
before were used. In this case, the real-time parameter 
estimation algorithm ran roughly 14 times faster than 
real time (1 second for a 14 second maneuver), because 
of the lower rate used for the recursive Fourier 
transform. 

In the final example, control surface inputs 
measured in flight during a piloted longitudinal/lateral 
tracking task were applied to a nonlinear F-16 
simulation 11 ' 14 with forward c.g. position (0.2 c). 

Figure 9 shows the control surface inputs. The 
maneuver was initiated from a steady trim condition at 
10 deg angle of attack, 10,000 ft altitude, and 
Mach 0.32. Simulated output data from the nonlinear 
simulation was corrupted with 20% white gaussian 
noise. The simulated noisy outputs are plotted in 
Figure 10. In this case, the lateral/directional linear 
model parameters were estimated in addition to the 
longitudinal model parameters from Eqs. (21)-(22). 


For the lateral/directional aircraft dynamics, the 
state vector x and input vector u in Eq. (1) are defined 
by 


x = [P p r <j>\ u = [S a 8,] (23) 

System matrices containing the model parameters are: 


A = 


N f 


sin a —cos a 


N„ 


B = 


0 

L s a 

N S a 

0 


L r 

N r 
tan 9 

Y *r 

L s, 

N Sr 

0 


— cos 9 


V, 


(24) 


(25) 


Figure 1 1 shows time histories of selected 
longitudinal and lateral/directional model parameter 
estimates, with dashed lines indicating the ±2<7 error 
bounds. Comparison values for the stability and control 
derivatives (plotted as solid lines) were calculated from 
the nonlinear simulation using central finite differences 
with a 1% perturbation size. Plots for the model 
parameters not shown were similar. As in the other 
examples, the performance of the real-time estimation 
algorithm was excellent, in spite of a relatively short 
(10 second) maneuver and low information content 
longitudinally. Parameter estimates converged to the 
finite difference values with appropriate standard error 
estimates. Cycle times for the short period and Dutch 
roll modes, based on the finite difference parameter 
values, were 3.8 seconds and 2.4 seconds, respectively. 
The real-time parameter estimation algorithm required 
4 seconds to converge to the finite difference values for 
the model parameters using the tracking inputs 
measured in flight and plotted in Figure 9. 

The sampling rate for the data was 80 Hz, and the 
recursive Fourier transform updates were done at 
40 Hz. For the same implementation and computer as 
before, the real-time estimation algorithm for the 
longitudinal and lateral/directional models together ran 
roughly 8 times faster than real time (1.2 seconds for a 
10 second maneuver). 
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Concluding Remarks 

A technique was developed for estimating linear 
model parameters in real time, using an equation error 
formulation in the frequency domain with a recursive 
Fourier transform. Simulation and flight test examples 
demonstrated that the method gives accurate real-time 
estimates of model parameters and standard errors. 

The examples also showed that the advantages of this 
approach include automatic noise filtering, robustness 
to high noise levels and data dropouts, fixed memory 
requirements regardless of the time record length, no 
bias parameters to estimate, and good performance for 
low information content in the data. All of these 
favorable characteristics, as well as low computational 
requirements, follow from analyzing the data in the 
frequency domain with a recursive Fourier 
transformation using fixed discrete frequencies within 
the frequency range for the dynamic motion of interest. 

The practical applicability of the method was 
demonstrated using a flight test data example, and a 
nonlinear simulation example using real flight test 
tracking inputs implemented by the pilot. Data 
information requirements for good parameter estimates 
were found to be low enough that they could be 
satisfied using ordinary pilot inputs measured in flight. 
The algorithm exhibited rapid convergence to accurate 
parameter values with standard errors that properly 
represented the accuracy of the parameter estimates. 
Typical convergence times were less than 1 cycle time 
of the dominant dynamic mode. No starting values 
were required for the parameter estimates, no tuning 
parameters had to be adjusted, and there was no 
temporal or spatial regularization during the parameter 
estimation. Parameter estimates and standard errors 
were based on measured data alone. The procedure 
was shown to have reasonable computational 
requirements, and ran much faster than real time, even 
when implemented in a high level language such as 
MATLAB. 

The technique could be used for dimensional or 
non-dimensional parameter estimation, and could also 
be used with general nonlinear models, as long as the 
model is linear in the parameters. All states and inputs 
must be measured, but this should not be a problem in 
modern aircraft with continuous automatic feedback 
control, for which the method is intended. 

The real-time parameter estimation technique 
studied in this work represents a fundamental building 
block for fulfilling the requirements of parameter 
estimation for adaptive or reconfigurable control. 
Future developments must focus on tracking rapid 


time-varying linear model parameters resulting from 
changes in the flight condition, stores, configuration, 
remaining fuel, or from various types of failures, wear, 
or damage. In addition, there must be some work done 
to address the question of excitation input design when 
insufficient information content in the data precludes 
accurate parameter estimates. 
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Table 1 F- 16 Linear Simulation Results, 


a () = 

1 deg, h 0 = 

10,000 ft, M 0 

= 0.37 

Parameter 

True 

20% Noise 

50% Noise 


V alue 

Estimate 

Estimate 



(Std. Error) 

(Std. Error) 

z a 

-0.600 

-0.624 

-0.602 



(0.047) 

(0.132) 

Z 'a 

0.950 

0.960 

0.986 



(0.016) 

(0.044) 

Z <L 

-0.115 

-0.104 

-0.134 



(0.017) 

(0.049) 

M a 

-4.300 

^4.197 

^1.021 



(0.136) 

(0.354) 

M q 

-1.200 

-1.238 

-1.200 



(0.045) 

(0.119) 

M 8 e 

-5.157 

-5.157 

-5.246 



(0.048) 

(0.130) 


Table 2 F18 HARV Flight Test Results, 
a G = 20 deg, h a = 24,100 ft, M C) = 0.34 



Batch 

Recursive 

Parameter 

Time Domain 

Frequency Domain 


Estimate 

Estimate 


(Std. Error) 

(Std. Error) 


-0.218 

-0.209 


(0.040) 

(0.084) 

z: 

1.047 

1.074 


(0.036) 

(0.052) 

Zz 

-0.057 

-0.041 

o s 

(0.046) 

(0.080) 

M a 

-0.649 

-0.509 


(0.044) 

(0.174) 

M 

-0.063 

-0.177 


(0.032) 

(0.107) 

M S, 

-1.257 

-1.415 

(0.063) 

(0.165) 
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Figure 1 Elevator Input 
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Figure 4 Simulated Measured Output, 50% Noise 
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Figure 2 Simulated Measured Outputs, 20% Noise 
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Figure 3 Parameter Estimation, 20% Noise 
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Figure 5 Parameter Estimation, 50% Noise 
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